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We address the balance of quantum correlations for continuous variable (CV) states. In 
T~H particular, we consider a class of feasible tripartite CV pure states and explicitly prove 

^ two Koashi- Winter-like conservation laws involving Gaussian entanglement of formation, 

C 2 ^ Gaussian quantum discord and sub-system Von Neumann entropies. We also address the 

class of tripartite CV mixed states resulting from the propagation in a noisy environment, 
CO and discuss how the previous equalities evolve into inequalities. 

o 
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1. Introduction 

Multipartite quantum correlations in continuous variable (CV) systems are valuable 

resources for quantum technology EEEEEIII Fully inseparable three- mode Gaussian 

states HEHni have been proposed to realize cloning at distanceEEH, and experimental 

schemes to generate multimode CV entangled states have been already suggested 
and demonstrated [i3|i4|i5|i6|i7|i8|i9:[ 

More recently, more general quantum correlations in bipartite Gaussian states 
have been analyzed: Gaussian quantum discord has been introduced and ex- 
perimentally investigated p2 | 23 | 24 ^ p^^, discrete variables several interesting rela- 
tions among the different measures of quantum correlations have been introduced 
|25|26|27|28|29|30|3i|32] ^j^j^ ^^^.^^ clarifying their different meaning and role 

The dimension of the Hilbert space does not appear to play a crucial role and, 
indeed, balance of correlations for CV systems has been derived in terms of the 
Renyi entropy In addition, balance of correlations in quantum measurements, 
and during the propagation in noisy channels, has been discussed both for discrete 
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and continuous variables |35l36|37|38|39|4o]^ 

Motivated by these results, in this paper we investigate whether some form of 
conservation laws for quantum correlations may be established for continuous vari- 
able systems and Gaussian measures of correlations. In particular, we consider a 
specific class of feasible pure CV tripartite states, and analyze in details the balance 
of Gaussian correlations among the bipartite subsystems. We found that Koashi- 
Winter-like relations, involving entanglement of formation, quantum discord and 
system entropies, may be generalized to continuous variables using Gaussian mea- 
sures of quantum correlations, at least for the class of feasible states that we took 
into account. The results are encouraging enough to suggest a direct experimen- 
tal verification and to foster investigations of correlation balance in more general 
Gaussian states. Our results may also taken as an argument in favor of the conjec- 
ture that Gaussian quantum discord is the quantum discord for Gaussian states 
and contribute to the ongoing discussion on the nature of quantum correlations in 
continuous variable systems^. 

The paper is structure as follows. In the next Section we introduce the class of 
feasible tripartite Gaussian states we are going to consider. In Section [3] we present 
the conservation laws for quantum correlations, whereas in Section |4] we consider 
the propagation in a noisy environment, and discuss how the previous equalities 
evolve into inequalities. Section [5] closes the paper with some concluding remarks. 

2. A class of feasible tripartite Gaussian states 

Let us consider three modes of a bosonic field, coupled through the interaction 
Hamiltonian 

H ^ jia'^ + 'f2b'^ c + h.c. . (1) 

The Hamiltonian H describes, e.g, two interlinked bilinear interactions taking place 
among three modes of the radiation field. In fact, it has been studied long ago 
for the description of parametric processes in x*^^^ media, assuming a suitable 
configuration satisfying phase-matching conditions. Interaction schemes described 
by H has been been realized using a single x*'^'' nonlinear crystal where the two 
interactions take place simultaneously f^^. The three-mode entanglement generated 
by H may be indeed exploited to realize optimal cloning at distance of coherent 
states lis!. Analogue interaction schemes may be realized in condensate systems USUI 
The Hamiltonian in Eq. ([T]) admits the constant of motion 

K{t) = NA{t) - Neit) - Ncit) = X(0), 

where Nj{t) denotes the average number of photons in the j-th mode. Starting from 
the vacuum |0) = |0)i(g)|0)2(g)|0)3 we have i^(0) = and thus TV^ W ^ NB{t)+Ncit) 
Vi. The expressions for Nj{t) can be obtained by the Heisenberg evolution of the 
field operators; upon introducing fl — v^|72p ~ |7iP we have 

Na = Nb+Nc, Nb^ '^''''7''' [cos - 1]2 , iVc = M! sin2(r!i) 
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The evolved state \T 
1 



Ut\0) reads as follows 



|T) = 



Si 



Nb 
+ Na 



n/2 



Nc 
1 + Na 



n\m\ 



\n + m,n,m) , (2) 



where Ut — exp {—iHt) is the evolution operator. The states |T) are Gaussian 
states, since they are generated from the vacuum by a bilinear Hamiltonian. Upon 
introducing the canonical operators qk and pk, k = A,B,C, and the vector R = 
{qA,PA,qB,PB,QC:Pc)^ 1 it is straightforward to prove that the mean values of the 
canonical operators are zero (R) = {T\R\T) = 0, whereas the covariance matrix 
(CM), whose elements are given by 

'^hk = ^(i^k, Rh}) - {Rk){Rh), 

can be written in the following block form: 

Da Oab Oac 
Oab Db Obc 



(3) 



where: 
Da 



Na + 



Df 



. Oac Obc 



Nb + 



Dc 



D( 



Nc + -]l (4a: 



Oab = V^b{Na + 1)^ Oac = ^/NcWI+1)F Ocb = V^bNc 1 , (4b) 



with Na = Nb + Nc 



Diag(l, 1) and P = Diag(l, —1). Since partial trace is a 



Gaussian operation the single-mode partial traces: 



QA = Ttbc 



\T){T\ 



and the two-mode ones 



qab = Trc 



|r)(r| 



qb = Ttac 



qac = Trs 



|T)(r| 



|T)(r| 



qbc = T^A 



\T){T\ 



\T){T\ 



are Gaussian states as well. The corresponding CM may be obtained by dropping 
the corresponding lines in Sy. The single-mode partial traces thus have diagonal 
CM CTk = Dk, corresponding to thermal states: 



k = A,B,C, 



^ (iVfc + l)P+i 

and the corresponding von Neumann entropies are given by 



(5) 



where: 



for a; > - . 
2 



(6) 
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The CM of the two-mode partial traces are given by 

( Da Oab\ ( Da Oac\ ( Db Obc\ 

and are already in the so called standard block form (to which any CM may be 
brought by means of local symplectic transformations) 

cr = A = Diag(a,a), B = Diag(5,fe), C = Diag(ci,C2) . 

For (Tab and (Tac we have ci — —ci whereas for ctbc we have c\ = C2. This means 
that QAB and qac are squeezed thermal states (STS) of the form: 

Nb 



QAB^S{T)v(Nc)®\G){Q\S\r), sinhV=— (8) 

QAC = S{T)v{NB)®\^){^\S\r), sinhV=-^^, (9) 

respectively, S(t) = exp{r(a^6^ — a6)} being the two-mode squeezing operator, 
whereas qbc corresponds to a mixed thermal state (MTS) that may equivalently 
expressed as (recall that N a = ^b + ^c)'- 

qbc^U{^)v{Na)®\^){^\U\cI,), cos2 0^— , (10) 

Nb + jVc 

QBC^U{<i,m{^\®v{NA)U\4>), cos^0^ ^1 , (11) 

where C/(0) = exp{0(a^6 — a6^)} is a bilinear mixing operator describing e.g. the 
action of a beam splitter. 

3. Balance of correlations 



In order to investigate the correlations between the modes of the state |T), we intro- 
duce the symplectic eigenvalues Aj_ of the two mode state Qh^ and the minimum 
symplectic eigenvalue }^^^^ of its partial transpose, which may be obtained asP^ 



y{hk) ^ 2-1/2 ^j^^ ^ 2/;,fe ± ^{h+Ik+2hk?~^Jhk , 

~^{hk) ^ 2-1/2 ^ + _ 2hk - Vih +h- llhk? - 4.Jhk , (12) 
respectively, where 

Jfc = det Dfc, Ihk = dot Ohk, Jhk = dct er^k ■ 

Note that: 

Sk = f (Vh) ■ (13) 



whereas the von Neumann entropies of the two-mode partial traces are given by 

Shk = fiX^t^) + fi\^-''^). (14) 
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The positivity of the state is equivalent to Aj_' > ^ and the positivity of the 
partially transposed state to A^'^'' > The bipartite Gaussian state ghk is thus 



entangled iff < aL''''^ < i In 

our case we have 



A 



(hk) 



- and A I 
2 + 



(hk) 



2 ^ 



(15) 



and thus the von Neumann entropy of any two-mode partial trace equals the von 
Neumann entropy of the remaining (single) mode, in formula: 



hk 



5, 



(16) 



where, in both Eqs. ( 15 ) and ( 16 1, we have h =^ k =^ j and h, k,j — A, B, C. 



For the symplectic eigenvalues of the partial transposes we have 

2 / , X 2 



:ab 



1 

71 



NA + \j +2Nb{Na + 1)+(nb + ^ 



(Na +Nb + l)^4NBiNA + l)+Nl; 



1/2 



X^(^ = X^^{Na^ Nb) , 



[BC 



(Nb - Ncf + Nb+Nc + ^-\Nb~ Nc\ ^{Na + 1) ^ - ANbNc 



(17) 
(18) 

1/2 

(19) 



from which it is easy to see that for any value of Nb and Nc we have 



< A 



{AB) 



< 



1 



< A 



(AC) 



< 



1 



A 



(BC) 



> 



1 

2 ' 



~(hk) 1 

i.e., QAB and qac are entangled states whereas qbc is separable. If < A_ < 5, 
i.e., Qhk is entangled, then the Gaussian entanglement of formation (EoF) is given 

by 



whereM 



Vhk = 



Ehkig) h{yhk) , 

iVh + v^KVhh - \hk\ + 1/4) - 2^\hk\Jhk 
iVh + VTk^ - Mhkl 



(20) 



with Jhk — det(cr/ife + |r2), and the symplectic form fl given by 

n = uj ® u) Lj 



1 

-1 

In our case, qab and qac are entangled and we have: 



1 



VAB = Vba = 



Nb 



2 l + Nc' 



VAC = ycA 



Nc 



1 



2 1 + Nb 



(21) 
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respectively. On the other hand, being the state qbc separable, yBc — Ucb ~ ^ 
and thus its EoF is zero. 

For MTS and STS the Gaussian quantum discords ^^{ghk) = T^hk and 
{Ohk) = T^kh (measurements performed on the first and second mode respec- 
tively) may be expressed as;^° 



T^hk = f [Vh) + f l^Vh - 
■Dkh^f(^/lk)+f(^/h- 



21/) 



hk\ 



2\Ihk\ 



l + 2^/I 



(22a) 
(22b) 



Starting from Eqs. (22 1 we can write a conservation law for the quantum corre- 



lations of the two-modes states Qhk- At first, we rewrite Eqs. (22 1 as follows: 

2\Ihk\ 



hk 



h X 



, (hk) 



= Sh + f Vh - 



■Dkh + h f aV^')) =Sk + f f - 



l + 2Vh 

2\hk\ 
l + 2JTk 



-/A 



Ahk) 



-/A 



ihk) 



(23a) 
(23b) 



where we used Eq. (|13|. Besides, for the state \2\, Eqs. (|15p say that / ( A 




{hk) 



and / ^A^'^-'^ = Sj, h ^ k ^ j, where we used the identity 
Furthermore, it is straightforward to verify that: 



A*;''-' = ./77 



1 + 2Jlk 



(24) 



with given in Eqs. (21 1. Summarizing, for the state \T) in Eq. ([2|, Eqs. (23 1 
leads to the following conservation laws 



'D/ik + Sj — Sh + Ekj 
T^hk + Shk — Sfi + Ekj 



(25a) 
(25b) 



with h k j and k, j = A, B, C. 

The above equalities are explicitly showing that Koashi- Winter-like conservation 
laws for quantum correlations may be written for continuous variable systems using 
Gaussian measures for quantum correlations, at least for the class of feasible states 
described by Eq. ^ . Eq. ( 25b I is the direct counterpart of the balance of correlations 
obtained for discrete variables, and thus we may expect, or conjecture, that it is 



of general validity for CV systems. On the other hand, Eq. (25a) is due to the 



specific properties of the state |T). Notice also that the validity of Eq. (25b) may 
represent an argument in favor of the conjecture that Gaussian quantum discord is 
the quantum discord for Gaussian states'^. Finally, we notice that both Eqs. (25a 



and ( 25b ) represent an independent check of the formula for the EoF in the case of 



non symmetric states'^. 
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4. Balance of correlations in the presence of noise 



In this section we investigate how, and to which extent, the conservation laws (251 
addressed in the previous section are modified by the propagation of the state in 
noisy channels. We assume that the three modes evolve through three identical 
uncorrelated noisy channels and that the Markovian approximation is valid. The 
propagation is thus described by the following Master equation 



7 

^*^2 

fc=l,2,3 



J2 {(Mh + i)/:M + iVth/:[4]} (26) 



where gt is the density matrix of the tripartite system described by the field oper- 
ators ai = a, a2 = b and = c, C[a]Qt = 2agta^ — a^agt — gta^a is the Lindblad 
superoperator, 7 is the overall damping rate, while A'th represents the effective num- 
ber of photons of the noisy channels Upon preparing the three modes in the initial 
state |r), the Gaussian nature is preserved during the evolution and the evolved 
CM St can be written asPl 

= e-^ St + (1 - e-") (27) 

where we introduced r — jt, Sx is the CM of the initial state given in Eq. ^ and 
5^00 = 5(1 + 2A'th)l4 is the asymptotic CM, I4 being the 4x4 identity matrix. 



By using the evolved CM ( 27 ) and its time-dependent local symplectic invariants. 



we can calculate the quantities appearing in Eqs. (23 1. In particular, now one has: 



/ (aL'''^) > and / (aV"'^) < S, (28) 



and, in turn, the equalities (251 changes into inequalities. The analytic expressions 
of the evolved conservation laws may be evaluated analytically, but they are quite 
clumsy and are not reported here explicitly. 



In order to study the evolution of the law ( |25a ) we define the following function: 



Aftfe (r) = [Vhk (r) + S, (r)] - [5^ (r ) + Ekj (r)] (29) 

where all the involved quantities are calculated starting from the results of the 
previous section but with the evolved state gt- As an example. Fig. [TJis a region- 
plot Ahk{T) as functions of r and S = Nb — Nc for given Nb and A^th^ depending 
on the values of the involved parameters, iS.hk{T) can be positive or negative. In 
particular, A^^ and IS.ac are always negative, A^a whereas d l^cA can change 
the sign (remarkably, this holds true also for A'th = 0). For the case of symmetric 
states, i.e., 5 = 0, one can write the following set of inequalities: 

Vab +Sc<Sa, Vba +Sc>Sb + Eac , (30a) 

T^Ac +Sb<Sa, -De a +Sb>Sc + Eab , (30b) 

Vbc +Sa>Sb+ Eca , VcB +Sa>Sc + Eba , (30c) 



where we used Ebc — Ecb ~ 0. 
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Fig. 1. (Color online) Balance of the correlations in the presence of dissipation and thermal noise. 
Region plots of A^j., h ^ k with h,k = A, B, C, as functions of 5 = Ng — Nq of t = ■yt. We set 
Nb = 1 and JVth = 0.2. Gray and white regions refer to the corresponding A^j, > and A^^ < 0, 



respectively. The red line represents the separability time Tsop: if r > r. 



scp then A^*' > | and 



the state Qhk{t) is no longer entangled (Ent.) and becomes separable (Sep.), i.e., Ej^i^ = 0. See the 
text for details. 



Analogously, to address the evolution of the conservation law ( 25b ) , we introduce 
the function: 



(KW) 
^hk 



(r) = [Vhk{T) + Shk{T)] - [Sh{T) + Ekj{T)] 



(31) 



As one can see in Fig. [ij where we plot Aj^^^-* (r) for the same choice of parameters 
as in Fig. [l] now one finds the following inequality holding for all the bipartitions, 
and any value of the interaction time, 



Vhk (r) + Shk {T)>Sh (r) + Ekj (r) , 



(32) 



h =/= k =/= j and h, k,j~A,B, C. Ineq. ( 32 1 generalizes to CV and Gaussian measures 
of correlation, the inequality discussed in '^^ for the discrete case. 



5. Conclusions 

In conclusion, we have proved that the balance of correlations originally investigated 
for three-qubit systems, involving entanglement of formation, quantum discord and 
single-system entropies is valid also for a feasible class of tripartite Gaussian CV 
states, upon using Gaussian measures of quantum correlations. Furthermore, in 
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Fig. 2. (Color online) Balance of the correlations in the presence of dissipation and thermal noise. 
Region plots of A^^^' , h k with h,k = A, B, C, as functions of <5 = Ng — Nq of t = 7t, for 
the same choice of the other involved parameters as in Fig. [T] Note that one always has A^^ > 
(gray regions). The red line still represents the separability time Tgcp. See the text for details. 



the presence of dissipation and thermal noise, the balance turns into inequalities 
between the previous quantities, depending on the actual values of the involved 
parameters. The results are encouraging enough to suggest a direct experimental 
verification and to foster investigations of correlation balance in more general Gaus- 
sian states. 
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